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ABSTRACT 

We conduct a series of comparisons between spectroscopic and photometric ob¬ 
servations of globular clusters and stellar models to examine their predictive power. 
Data from medium-to-high resolution spectroscopic surveys of lithium allow us to in¬ 
vestigate first dredge-up and extra mixing in two clusters well separated in metallicity. 
Abundances at first dredge-up are satisfactorily reproduced but there is preliminary 
evidence to suggest that the models overestimate the luminosity at which the surface 
composition first changes in the lowest-metallicity system. Our models also begin extra 
mixing at luminosities that are too high, demonstrating a significant discrepancy with 
observations at low metallicity. We model the abundance changes during extra mixing 
as a thermohaline process and determine that the usual diffusive form of this mech¬ 
anism cannot simultaneously reproduce both the carbon and lithium observations. 
Hubble Space Telescope photometry provides turnoff and bump magnitudes in a large 
number of globular clusters and offers the opportunity to better test stellar modelling 
as function of metallicity. We directly compare the predicted main-sequence turn-off 
and bump magnitudes as well as the distance-independent parameter AMy 
We require 15 Gyr isochrones to match the main-sequence turn-off magnitude in some 
clusters and cannot match the bump in low-metallicity systems. Changes to the dis¬ 
tance modulus, metallicity scale and bolometric corrections may impact on the direct 
comparisons but AMy which is also underestimated from the models, can only 

be improved through changes to the input physics. Overshooting at the base of the 
convective envelope with an efficiency that is metallicity dependent is required to 
reproduce the empirically determined value of AMy 

Key words: stars:abundances, evolution, interiors. Population II. 


1 INTRODUCTION 

Globular clusters (GCs) host mono-met allicQ, almost coeval 
stellar populations that have been studied extensively, both 


* E-mail: angelou@mps.mpg.de 
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Johnson & Pilachowskil l2010h. M22 jMarino et al 

l2009f). M54 

Carretta et al.ll2010ah. M2 ilVong et al.||2014r) 


photometrically and spectroscopically. In spite of their mu l- 
tiple populations fsee lGratton. Garretta fc Bragagliall2012l '). 
these systems offer well constrained tests for stellar evolu¬ 
tion theory. Their colour magnitude diagrams (GMD, hence 
photometry) provide two clear evolutionary diagnostics for 
the stellar models - the main-sequence turn-off (MSTO) 
magnitude and the magnitude of the bump in the red gi¬ 
ant branch (RGB) luminosity function (LF bump or bump 
hereinafter). The MSTO magnitude serves as the classic in- 
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dicator for cluster age, and by fitting this region of the CMD, 
modellers can take solace in the fact their codes provide a 
reasonable approximation to the stellar physics during the 
early phases of evolution. The bump magnitude, on the other 
hand, reveals information about the depth of first dredge-up 
(FDU) and stellar mixing - it is a probe of internal processes 
that are significant in later stages of evolution. Spectroscopy 
complements (and vice-versa) these photometric studies by 
providing further quantitative details of the internal mix¬ 
ing processes. It is our interest in stellar abundances and 
mixing processes, particularly during the RGB ascent, that 
motivate this study. 

Mixing in RGB stars is characterised by two distinct 
events; the well-understood FDU and an additional mixing 
episode that operates over a longer timescale during a more 
advanced phase of RGB evolution (“extra mixing” here¬ 
inafter). Only mixing dur ing FDU is predicted by standard 
stellar theory (||be3[l963). In low-mass stars (M < 1 Mq), 
deep convective motions develop once the star becomes a gi¬ 
ant and these can penetrate into regions that have previously 
experienced partial hydrogen burning. Material enriched in 
'*He, and is mixed through the convection zone in¬ 
creasing the prevalence of these nuclei at the stellar surface. 
Conversely, ^Li and abundances decrease as they are di¬ 
luted into the extending envelope. In solar-metallicity stars 
the ratio falls from the solar value of « 90 to « 30 

after the FDU, whilst ’^Li is depleted by a factor of « 20. 
Up to this point the theoretically predicted changes are in 
good agreement with observations. 

It is clear from observations that after FDU the 
surface co mposition is fu r ther a ltered during RGB 
evolution jGilrov fc Brown _ 1991; iGratton et al. _ 200^ 


Snhth&_^artell gOOj;^ Shetrond 20031: 1 Weiss fc Charbonn3 


20041 : Martell. Smith fc Brilev 


This second mixing 
event sees the cycling of hydrogen burning products into 
the convective envelope. From a post-FDU value of approx¬ 
imately 30, the ratio reduces to « 15 in solar- 

metallicity stars and to the equilibrium value of « 4 in 
metal-poor stars. Lithium is essentially destroyed at all 
metallicitie^. These changes are not predicted by standard 
stellar theory. 

The onset of extra mixing seemingly coincides with 
the bump in the luminosity function of GCs. The mixing 
is therefore associated with the advance of the hydrogen¬ 
burning shell into the composition discontinuity left be¬ 
hind by the deepest extent of FDU. The internal pro- 
cess(es) responsible for the surface changes are thought 
not to manifest themselves until after the bump because 
the composition discontinuity, and hence discontinuity in 


rier to any extra mixing process ( 

Sweigart & Mengelll 19791: 

ICharbonnel. Brown & Wallerstein 

Il998h. 


In this study we are concerned with how well stellar 
models reproduce the key properties (magnitude/luminosity 
onset, abundance changes) of several evolutionary features 
experienced during the late main-sequence and RGB phase 
of evolution. These include: 


the MSTO magnitude. 


^ Save for the small number of Li-rich giants discovered. See 
dP’Orazi et al.l^2015a^ for example. 


• the LF bump magnitude, 

• the difference between the MSTO and bump magni¬ 
tude, 

• the FDU magnitude, 

• the abundance changes associated with FDU, 

• the abundance changes associated with extra mixing. 

Observations of lithium in GC stars are used to 
constrain FDU and the extra mixing event. Match¬ 
ing the magnitude of the LF bump has, in the past, 
been used as a measure of the adequacy o f stel- 


lar 

models (e.g., iKing. Da Gosta & Demarouel 19851: 

Fusi Pecci et al.lll990l: lAlongi et al.lll99ll: ICassisi & Salarig 

1997 

: Gassisi. degl’Innocenti & Safaris! 1997: Zoccali et al. 

199C 

: ISalaris. Gassisi & WeissI l2002l: iRiello et al.l 

2003|: 

Biork & Chaboved l2006l: iMeissner & WeissI 

20061: 

Michaud. Richer & RichardI l2010l: Di Cecco et al.l 

20ld; 

Gassisi et al.l |2011^. Here we draw upon some of the 


methodology of these studies and compare to recently 
available large data sets. 

We first turn our attention to spectroscopic surveys of 
lithium abundances in GCs and focus, in particular, on two 
clusters that are separated in metallicity by factor of ten. 
These systems provide a detailed examination of FDU and 
extra mixing. In order to better sample the metallicity dis¬ 
tribution, we utilise Hubble Space Telescope (HST) pho¬ 
tometry to determine how well stellar models reproduce the 
MSTO magnitude, bump magnitude and the difference be¬ 
tween the two (which is independent of distance and red¬ 
dening) in a large number of systems. These three compar¬ 
isons allow us to investigate the uncertainties introduced 
from the metallicity scale, distance modulus and bolometric 
corrections. We determine how changes to the input physics 
impact upon the level of agreement between observations 
and models. We conclude with a spectroscopic analysis of 
extra mixing and the ability of the thermohaline mecha¬ 
nism to simultaneously account for the depletion of carbon 
and lithium as a function of luminosity. Such a demand is a 
stringent test of any extra mixing process. 


2 OVERVIEW OF THE STELLAR MODELS 

The stellar models in this work are calculated with MON¬ 
STAR (t he Monash version of the M t. Stromlo evolution 
code; see ICampbell fc LattanzT^ l2008f l. Opacities are pro¬ 
vided at the high te mperature end by the O PAL Rosseland 
mean opacity tables dlglesias fc Rogerall996l L Our standard 
procedure is to employ tables that are variable in G and 
O content (OPAL type-2 ) based on the so l ar he avy ele¬ 
ment mixture specified by Grevesse fc Noei3 (Il993l . GN93). 
The code also utilises additional (OPAL type-1) tables with 
a fixed metal distribution when required. For this work 
we have gener a ted fix ed metal distribution tables with the 
lAsplund et al.l ll2009l . AGSS09) solar mixture, and alpha- 
element enhancements of [Q/Fe]= 0.2 and [a/Fe]= 0.4 of 
the AGSS09 abundances . Below 10"^K opacity tables from 
iLederer fc Aringe^ 200^ with variable C and N content are 
used fsee lCampbel l2007l and iGampbell fc Lattanzid l2008l 
for further details). 

MONSTAR ordinarily only follows those species that 
are significant energetically. A seven species network (^H, 
®He, "^He, ^^N, ^®0, as well as a seventh pseudo-element 
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Reaction Rate 

Source 


^H(p, e+i/e)^H 

Harris et al. flQSS) 


3He(3He, 2p)4He 

Cauffhlan Sz Fowler i 

(1988) 

3He(4He, 7)’'Be 

Cauffhlan Sz Fowler I 

(1988) 

I2c(p, 7)13n 

Cauffhlan Fowler I 

F 

i^N^p, 7)150 

Champagne (2004, private comm) 

^Be(e“, 

Reaclib electron capture database 

7Be(p, 7)8B 

Aneulo et al. fl999^ 


l’Li(p, 4He)'iHe 

Descouvemont et al. 

(2004) 

i3C(p, 7)1'‘N 

Aneulo et al. ('1999^ 



Table 1. Key reaction rates used in MONSTAR. 


The equation for the diffusion coefficient is: 

where ip = (9 In p/9 In/i)p,T, 5 = —(91np/91nr)p,;i, 
= (91np/91nP), Vad = (91nr/91nP)ad, V = 
(91nr/91nP), K is the thermal diffusivity and Ct is a di¬ 
mensionless free parameter. In this diffusive theory, Ct is 
related to the aspect ratio, a, of the thermohaline fingers 
(assumed to be cylindrical) by 


that ensures baryon conservation) is sufficient to include 
feedback on the structure from the nuclear energy gener¬ 
ation. We have extended the network for this study so that 
the evolution of ^Be, ^Li and are now followed; with 
the necessary changes to temporal and s patial resolution 
criter ia to follow these fragile species (see iLattanzio et al.l 
I 2 OI 5 I for the importance of these criteria). The key reac¬ 
tions and source of each adopted rate can be found in Ta¬ 
ble [T] Note that we employ the ^'^N(p, 7)^®0 reaction rate 
provided by Champagne (priv ate communication) wh ich is 
consistent with that givenbyjAdelbergeret_^L| (1201 il l and 
the LUNA collaboration, [^almerini et al. 1 201llf discuss the 
consequences of the new rate for low-mass stellar evolution, 
including extra mixing. 

Convective energy transport in MONS TAR is treated 
accor ding to the mixing length theory (MLT, iBohm-Vitensd 
Il958lf . Mixing of th e chemical species is calculat ed using a 
diffusion equation (ICamobell fc Lattanzid [2008l L With an 
assumed helium content of Y= 0.2485 and heavy element 
mixture specihed by GN93, the MLT parameter ompt = 
1.75 best reproduces the solar model at the Sun’s current 
age. A value of omlt = 1.69 is required when the AGSS09 
abundances are adopted. In the models presented here, we 
employ the Schwarzschild criterion and where specified, al¬ 
low for non-locality by including a prescript ion for diffu¬ 
sive o vershoot. We follow the procedure of iHerwig et al.l 
lll997tl wh o parameterised the numerical sim ulations of con¬ 
vection bv iFrevtag. Ludwig fc SteffenI (Il996l i. Those results 
indicated an exponential decay in velocity of the overshoot¬ 
ing material. In analogy to the pressure scale height, Hp, a 
‘velocity scale height’, Hv, is defined such that 


= fosHp 


( 1 ) 


where fos is a scaling factor which we vary between fos = 
0.0—0.1 in this study. The resulting equation for the diffusion 
coefficient is then 


Do3 = Do e 


( 2 ) 


where Dq is the diffusion coefficient at the last convective 
point and z is distance from the convective boundary. 

Our implementation of the ther mohalin e mech anism 
uses the formulation developed by [Uln^ (Il972l i and 
iKippenhahn. Ruschenplatt. fc Thomas! ( 1980l i. in which 
thermohaline mixing is modelled as a diffusive pro¬ 
cess. Thi s prescription has b e en em ploy ed in previous 
work by Gharbonnel fc Zalmj_j[ 20n7bj|al^ Stancliffe et al.l 
ll2n09h. IS tancliffj j2010h . rCharbonnd fc Lagardel l|2010ll . 
and Angelou et al.l J 2 OI 1 I . l2012fl 


822 
Gt = -TT a . 


(4) 


The mechanism is elegant in that the depth of mixing 
is set by the stellar structure resulting in only one free 
parameter. An empirically derived value of Ct = 1000 
can reproduc e abundance patte r ns in both globular 
cluster stars j[AjrgglouetjilJ [ 20111 . l2012fl and field stars 
llCharbonnel fc Lagardd l2010h as well as the dichotomy 
between car bon-normal and ca rb on-enhanced metal - 
poor stars dsta ncliffe et al.l l2009l i. iDenissenkovI (l2010l ~) 
and ICantiello fc Langerj 1 201^ ^efer a lower value of 
Ct = 12 on theoretical grounds, which is also supported 
by 3D hydrodynamical models that suggest the mixing 

is _ inefficient on the RGB dPenissenkov fc Merrvfieldl 

201l| _ [TraxJei\_Geraud_j^tellmaclTl 1201 ll : 


Brown. Garaud fc Stellmach 2013l i. Parametrising ex¬ 


tra mixing in this form and comparing to observations 
can still tell us much about the transport of material 
in the stars. Similar conclusions could be reached had 
we elected to include a ph e nome nological mixing model 
dPenissenkov fc VandenBerd l2003lL Note that in such 
models the depth of mixing is usually specified by a 
constant shift in either mass or temperature from the 
hydrogen-burning shell. An additional free parameter sets 
the mixing speed. This is a less physically motivated but 
similar configuration to that used here. 

In order to compare to observations, we convert our 
stellar models from luminosity space to absolute visual mag¬ 
nitude [Mv) throughout. This requires a V band bolomet- 
ric correction for each st ellar model calculated. Theoretical 
model atmospheres from ICastelli. Gratton fc Kurucj dl997l . 
ATLAS9) provide the necessary tables of bolometric correc¬ 
tions and we determine the most appropriate value by using 
a cubic spline to interpolate in composition, surface gravity 
and Teff. 


3 FIRST ANALYSIS OF THE 
SPECTROSCOPIC DATA 

3.1 Observations of Lithium in Globular Clusters 

The lithium abundances used to constrain our models come 
from medium-to-high resolution surveys of GCs. Although 
far more complex than their classic simple stellar popula¬ 
tion archetype, GCs are still useful testbeds of stellar the¬ 
ory due to their well populated colour-magnitude diagrams 
(CMDs). These systems host multiple stellar populations 
spanning relatively small differences in age and are char¬ 
acterised by internal variations in their light-element (e.g., 
C, N, Na, O, Al, see iGratton. Carretta fc Bragaglial [ 2 OI 2 I 
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and references therein for a recent review on multiple pop¬ 
ulations) and ^He content (iGratton et al.lBoid . l201h . I2OI2I : 
iMilone et al.l[20l3 i. 

In addition to these primordial abundance varia¬ 
tions, the constituent stars undergo in situ composi¬ 
tion changes such as those experienced on the RGB. In 
situ mixing manifests itself as a function of luminosity 
and is easily discerned. The carbon and nitrogen abun¬ 
dances, in particular, have been extensively measured. The 
I2c/i3c j-atio is also a useful tracer of mixing. It scru¬ 
tinises FPU robustly JPe arborn. Bolton fc EggletonI 1 19751 : 
iTomkin. Luck fc Lambertlll976l : ICharbonnellll994h . but sat¬ 

urates rapidly once extra mixing begins. Measuring 
requires spectra with medium-to-high resolution and signal- 
to-noise ratios. 

Lithium is a very useful probe of mixing and a sensitive 
gauge of temperature because it is destroyed at 2 MK. It 
can therefore provide an indication of the mixing efficiency 
of both FDU and the extra-mixing event. It also has the 
advantage that its abundance determination is reliable: it 
is derived from the the Li l resonance doublet at 6707.78A, 
rather than from molecular bands which is the case for G 
and N. 

We have compiled data from four medium-to-high res¬ 
olution studies that focus on lithium abundances in glob¬ 
ular cluster giants. Data from all surveys have typical un¬ 
certainties of roughly ±0.1 in A(Li) and ±0.03 mag. The 
aforementioned internal (i.e., star-to-star) errors in Li abun¬ 
dances are due to a combination of uncertainties in equiv¬ 
alent width measurements, continuum placement, signal-to- 
noise ratios of the spectra and atmospheric parameters (ef¬ 
fective temperatures, micorturbulence, gravity and metallic- 
it y^ with the first contr ibution being the largest). We refer 
to (p’Orazi et al.l (l2014l i for an extensive discussion of this 
topic. Our selected sample of clusters spans a large range 
in metallicity and we expect to identify any such trends 
present. In Table [2] we list the clusters with Li abundances 
determined and t heir general prope rties. Th e meta ll icities 
are as provided in lNataf et al.l il2013f l and the lHar^ (Il996l l 
catalogue (2010 edition). 

The most metal-rich cluster in our sample, M4 ([Fe/H] 
= —1.10), and the most metal-poor, NGC 6397 ([Fe/H] 
= —2.10), have subgiant bran ch lithium abundan ces con¬ 
sistent with the Spite plateau (ISoite fc SpitdfT982l . see also 
Figure 1); a point discussed by the authors of the original 
surveys. In these two clusters the first instance of lithium de¬ 
pletion marks the onset of FDU. The difference in metallicity 
causes stars in M4 to begin FDU approximately AMv ~ 0.3 
magnitudes fainter than those in NGC 6397. This is because 
in metal-rich stars, the additional low-ionisation metals (i.e., 
Ca, Na, K, and Al) contribute to the higher opacity. All 
other things being equal, the metal-rich stars: 

(i) begin the inward migration of the convective envelope 
at a lower luminosity; and 

(ii) develop deeper convective envelopes. 

compared with metal-poor stars. As one might expect from 
deeper dredge-up, stars in M4 exhibit a post-FDU Li abun¬ 
dance that is generally lower than those in NGC 6397. How¬ 
ever it is not certain that we are detecting a metallicity ef¬ 
fect because systematic uncertainties may exist between the 
two studies. The upper range of the Li abundances in M4 
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Figure 1. Top Panel: A(^Li) as a function of absolute magnitude 
{My) for the clusters specified in the legend. Bottom Panel: As 
above but with an offset applied to the magnitude so that the LF 
bumps coincide. In this panel we omit data from NGC 2808 as 
none of the stars were observed in the luminosity range where Li 
depletion via extra mixing is expected to occur. 


is consistent with the stars in NGC 6397. Furthermore, the 
significant change in abundances during FDU is due to the 
transition from shallow surface convection to a deep convec¬ 
tive envelope. A slightly deeper convective envelope may not 
necessarily translate to a detectable difference in the surface 
abundance. The two studies found the same pre-FDU abun¬ 
dances which might suggest that the systematic errors are 
not signihcant and the differences in the mean pre-bump Li 
abundances are indeed due to metallicity. Systematic uncer¬ 
tainties that exist between surveys, however, can only truly 
be minimised through a homogeneous study of both clus¬ 
ters (i.e., identical instrument, line lists, codes, methodol- 
og y). Interestingl y ^ thos e clu sters self-consis t ently a nalysed 
by iD’Orazi et al.l i2014h and iD’Orazi et al.l (l2015bll Span a 
smaller metallicity range (see Table El) but share the same 
mean pre-bump Li abundance once observational uncertain¬ 
ties are taken into account. 

The two clusters show noticeably different gradients of 
lithium depletion during FDU. Stars in M4 complete FDU 
over a much smaller luminosity range than those in NGC 
6397. The abundance trend reflects the rate of advance of 
FDU which is different because the greater opacity allows 
the convective envelope to penetrate faster in the metal-rich 
regime. 

Figure [T] also provides insight into the extra mixing pro¬ 
cess. The GCs in our sample have been surveyed across 
a large luminosity range that includes the LF bump. Be¬ 
cause the depth of FDU is metallicity dependent, so too 
is this secondary mixing event. Metal-rich stars begin ex¬ 
tra mixing at fainter magnitudes because the hydrogen shell 
is not required to advance as far before it encounters the 
homogenised region and removes the ^ inversion that in- 
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ID 

Messier ID 

Relative 

Age 

Luminosity 

{Mv) 

[Fe/H] 

RGBi^un^p 

{My) 

FDU 

{Mv) 

Number of 
Targets 

Source 

NGC 6121 

M4 

0.91 

-7.19 

-1.10 

0.45 

3.65 

87 

Mucciarelli et al. 12011') 

NGC 2808 

- 

0.74 

-9.39 

-1.18 

0.65 

- 

68 

D’Orazi et al. f2015b'l 

NGC 5904 

M5 

0.81 

-8.81 

-1.29 

0.50 

- 

99 

D’Orazi et al. ('2014') 

NGC 362 

- 

0.74 

-8.43 

-1.30 

0.57 

- 

67 

D’Orazi et al. ('2015b') 

NGC 6218 

M12 

0.92 

-7.31 

-1.37 

0.78 

- 

63 

D’Orazi et al. (2014) 

NGC 1904 

M79 

0.87 

-7.86 

-1.58 

0.29 

- 

47 

D’Orazi et al. f2015b'l 

NGC 6397 

- 

0.99 

-6.64 

-2.10 

0.16 

3.3 

454 

Lind et al. f2009^ 


Table 2. General pro perties of the clusters from which spectroscopic data is utilised in this work. Relative ages are as defined in 
IPe Angeli et all ll2nnfil i. 


hibits the mixing process. We note that such a metallicity 
trend is not evident from Figure [T] This clear in Table [2] 
which shows the statist ically determined LF bump magni¬ 
tude (|Nataf et al.ll201^ ') for each cluster sorted by metal¬ 
licity. The metal-rich clusters do not necessarily begin ex¬ 
tra mixing at fainter magnitudes as theory predicts. Uncer¬ 
tainty in the metallicity □ and distance modulus determina¬ 
tions and the role of multiple populations each contribute 
to the observed behaviour. Such disagreement with direct 
comparison of the evolutionary events is a common theme 
throughout this paper. 

In the bottom panel of Figure [T] we have applied a mag¬ 
nitude offset to each cluster so that the magnitude at which 
extra mixing begins is common to all clusters. In all clusters, 
irrespective of metallicity, extra mixing depletes lithium at a 
similar rate. Agreement is further improved when one con¬ 
siders that the abundances of the brightest three stars in 
NGC 6397 are upper limits. Such strong agreement is per¬ 
haps not unexpected given the fragility of lithium, but it also 
highlights a clear property (and constraint) of the physics 
that drives the mixing during this epoch. 


3.2 Surface Composition Changes During the Red 
Giant Branch Mixing Events 

3.2.1 Sheddinq New Liqht on Previous Work with fC/Fel 
and [N/Fe] 

When tr ying to match the extr a mixing event in M92 ([Fe/H] 
= — 2.21. lAngelou et al.l (l2012l l found that their models un¬ 
derestimated the magnitude (hence overestimated the lumi¬ 
nosity) of the LF bump by AMy ~ 0.7 mag. They were 
required to artificially deepen FDU significantly 0 in their 
calculations to match the photometrically determined mag¬ 
nitude of the LF bump. An extension of the convective enve¬ 
lope was also required to match the LF bump in M15 (also 
[Fe/H] = —2.2). In these clusters, it is unclear if the onset 
of extra mixing coincides with the photometric bump. One 
interpretation of the data is that surface abundance changes 
begin before the LF bump. Possible reasons for the discrep¬ 
ancy were given as uncertainties in the spectroscopy (e.g., 
combining data sets, determining abundances from molecu¬ 
lar bands), difficulties in determining the luminosity of the 


PT 


3.0 
2.5 
2.0 
.1.5 r 


1.0 

0.5 

0.0 


k A^ AA 



: A M92: [C/Fe] 


Lli 


NGC 6397: A(^Li) 

I .... I .... I .... I 


l‘ 

I 

li 


2 1 

Mv 


0 -1 -2 


0.5 

0.0 


-0.5 


- 1.0 


-1.5 


- 2.0 


3.0 

2.5 

2.0 

Sl.5 

< 

1.0 

0.5 

0.0 


0.5 
0.0 
-0.5 
- 1.0 
-1.5 
- 2.0 

5 4 3 2 1 0 -1 -2 


_ I I M M I M M I M M I M M I M M| I M M I M M I M M_ 



: A NGC 5466: [C/Fe] 

— , • 

i • NGC 6397: A(7Li) ' 


I I I I I I I I I I I I I I I I I I I I I I I I I il I I I I I 1*1 I I I I I I I i' 


Mv 


Figure 2. Left Axis (both panels): Black circles denote A(^Li) as 
a function of magnitude for stars in NGC 6397 ([Fe/H]= —2.1). 
The dashed vertical line corresponds to the photometrically de¬ 
termined LF bump of the cluster (My = 0.163). Right Axis (top 
panel): Lavender triangles denote [C/Fe] as function of magni¬ 
tude for stars in M92 ([Fe/H]= —2.2). The dotted vertical line 
corresponds to the photometrically determined LF bump of the 
cluster (Mv = 0.016). Right Axis (bottom panel): Lavender trian¬ 
gles denote [C/Fe] abundance as function of magnitude for stars 
in NGC 5466 ([Fe/H]= —2.2). The dotted vertical line corre¬ 
sponds to the photometrically determined LF bump of the cluster 
( My = —0.075 ) . The LF bump for each cluster was determined 
bv iNataf et al.l ll2013lV 


^ Whilst there may be some uncertainty in the metallicity deter¬ 
mination of these systems, it is a robust result that NGC 6397 and 
M4 differ by approximately a factor of 10 in their metal content. 
^ The depth of FDU was extended from m= 0.368 Mq to m = 
0.320 Mq, where m is the enclosed mass. 


bump at low metallicity, or that extra mixing had initiated 
before the LF bump (which would prove to be a serious issue 
for stellar evolution). The authors noted that homogeneous 
lithium data would reveal the true behaviour of the cluster. 

The three low-metallicity clusters studied by Angelou 
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et al. (2012) [M92, M15 and NGC 5466, [Fe/H] = -2.2] are 
yet to have their lithium abundances systematically mea¬ 
sured. They do, however, have comparable metallicity to 
the surveyed NGC 6397 ([Fe/H] = —2.1). In the upper panel 
of Figure [2] we plot [C/Fe] from M92 (right axis, lavender 
triangles) and A(Li) from NGC 6397 (left axis, black cir¬ 
cles) as functions of absolute visual magnit ude. Data for 
M92 are taken from ISrnith fc Martell J200d) who applied 
offsets to t he studies by Carbon et al.l (|l982l '). iLanger et al.l 
lll986tl and iBellman et al.l l|200lh in order to remove sys¬ 
tematic differences in abundance scales. Combining results 
from different studies, together with inferring abundances 
from molecular bands, may lead to uncertainties of up to 
0.3 dex in the [C /Fe] data. Data for NGC 6397 are from 
iLind et "^ ll200 9lf and include their re analysis of the sur- 
vev bv lconzalez Hernandez et aP (l2009l ). The dotted laven¬ 
der line indic ates the magnitud e of the LF bump in M92 
{Mv ~ 0.016; iNataf et al.l 201^ 1 and the black dashed line 
the magnitu de of the LF bump in NGC 6397 {Mv = 0.163; 
iNataf et al.l 12011^ . The magnitude at which extra mixing 
begins in the massive GC M92 is unclear from [C/Fe] from 
spectroscopy. The lithium decrease in NGC 6397, on the 
other hand, has a well defined starting magnitude that cor¬ 
responds to its photometrically determined bump and the 
bump of the similarly metal-poor M92. In the lower panel 
we compare data from the less massive cluster NGC 5466. 
As per the panel above, [C/Fe] is denoted by l avender tri¬ 
angles with the scale provided on the right axis. INataf et al.) 
(l2013tl determined the magnitude of the LF bump in this 
cluster to be My ~ —0.075. In NGC 5466 and NGC 6397 
the respective magnitudes of the LF bump and onset of sur¬ 
face abundance changes (due to extra mixing) agree quite 
well. 


Unlike M92, for which the data are a combination of 
multiple surveys, data for NGC 54 66 are homogeneous (t hey 
are taken from a single study bv IShetrone et aDl2010h . It 
may be that the combination of multiple data sets clouds 
the true cluster behaviour. If the uncertainty surrounding 
M92 is purely due to observational spread, then results from 
homogeneous surveys such as APOGEE will shed light on 
this issue. However, we note that both NGC 6397 and NGC 
5466 ar e mu ch less massiv e than M92 and I Garrett a et al.l 
ll2010bl l and iD’Orazi et al.l ll2014h have demonstrated how 
clusters with greater mass can exhibit a greater spread in 
their primordial abundances. 


Because the Li abundances suggest that the beginning 
of extra mixing in low-metallicity clusters does indeed co¬ 
incide with the magnitude of the photometrically derived 
LF bump, then it i s clear that the models presented in 
lAngelou et al.l 1 I 2 OI 2 II were underestimating the magnitude 
of the bump and of the onset of extra mixing. This is irre¬ 
spective of the chosen extra mixing mechanism and depen¬ 
dent on the physics of the stellar codes. Such an inconsis ¬ 
tenc y was identified b y King. Da Costa fc Demaraud (Il985l 'l 
and iFusi Pecci et al.l l|l990l L We explore the extent of this 
discrepancy in §4. 



NGC 6397 

M4 

Mass {Mq) 

0.80 

0.812 

z 

0.00011 

0.0011 

X('‘He) 

0.24989 

0.2489 

X(iH) 

0.75 

0.75 

X(12c) 

1.99 X 10“® 

1.99 X 10-"! 

X(14n) 

5.83 X 10-® 

5.83 X 10-® 

X(160) 

4.82 X 10“® 

4.82 X 10-1 

X("Li) 

9.39 X 10-1° 

9.39 X 10-1° 

«MLT 

1.75 

1.75 

Ct 

1000 

1000 

TO Age (Gyr) 

12.0 

12.0 

Opacity Mixture 

AGSS09 

AGSS09 


Table 3. Model details used to fit the clusters NGC 6397 and M4 
in Figure[3] Note that we require the same initial X(^Li) despite 
the factor of ten difference in metallicity. 


3.2.2 Comparison with Models: Elucidating the Red Giant 
Branch Mixing Events at Globular Cluster 
Metallicities 


Our sample of GCs includes the metal-poor cluster, NGC 
6397 ([Fe/H] « —2.10), and the metal-rich, NGC 6121 (M4, 
[Fe/H] « —1.10), which have been surveyed across a lumi¬ 
nosity range that covers FDU through to the extra mixing 
event. Li abundance determinations for NGC 6397 are pre¬ 
sented in the top row of Figure [3] and for M4 in the bottom 
row. The data are plotted both as a function of absolute 
visual magnitude (left panels) and as a function of indepen¬ 
dently determined luminosity (right panels). The two bright¬ 
ness systems are employed as a check on possible systematic 
errors in our method (see below). We calculate stellar mod¬ 
els for each cluster (lavender curves) with the parameters 
specified in Table [3] as a first approximation. 

The models for NGC 6397 underestimate the magni¬ 
tude of both mixing events. The difference in magnitude 
at FDU (« 0.3 mag) is not as pronounced as for the LF 
bump (« 0.7 mag). In the models for the metal-rich clus¬ 
ter, M4, the onset of FDU is consistent with the observations 
whilst a discrepancy is present at the LF bump (« 0.4 mag). 
It is well documented that theoreti cal models underesti¬ 
mate the magnitude of the LF bump (iFusi Pecci et al.lll990l: 


Biork fc Chaboverll2006l : IPi Cecco et al.l 2010l : ICassisi et al.l 
2 OIIII . however this is the first time a disagreement at FDU 


has been identified. The recent availability of lithium obser¬ 
vations, which probe both RGB mixing events, allows for 
an investigation of FDU across a range of metallicities. Pre¬ 
vious investigations of mixing during the RGB have relied 
on [C/Fe] and [N/Fe] which are of limited use for the FDU 
event because these species change very little during FDU at 
metallicities typical of GCs. Furthermore, their abundance 
determination (usually from molecular bands) is less robust 
than for Li. We have seen that in the lowest metallicity 
clusters there is an intrinsically large spread in C and N 
which m akes identifying the onset of mixing by spectroscopy 
difficult. ICharbonnell (Il994l . Il995l l has used both 


and Li to test extra mixing due to rotational instabilities. 
Model comparisons to M4 and Halo stars with NGC 6397 
metallicity were presented with their predicted FDU mag- 
nit udns_com£ajableor slightly brighter than ours. Models 
bv IPlacco et al.l ll2014lf determine that reproducing the on- 
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Figure 3. Top panels: A(Li) as a functi on of absolute ma gnitu de (left) and independently deter mined luminosity (right) in the globular 
cluster NGC 6397. Data are taken from iLind et al.l ll2009l l and I Gonzalez Hernandez et al.l ||2009| V Bottom panels: A(Li) as a function of 
abso lute magnitude (left) an d independently determined luminosity (right) in the globular cluster M4. Data for this cluster are taken 
from iMucciarelli et al.l ll201lh . In each panel the solid lavender curve denotes a first approximation model for the respective cluster (model 
parameters are listed in Table [3]) • 


set of extra mixing is also a problem for carbon-enhanced 
metal-poor stars. 


3.3 Bolometric Corrections 

It is concerning how much the models overestimate the 
brightness of the mixing events, especially the LF bump 
at low metallicity. It is prudent to first ascertain whether 
the magnitude offset between the theoretical and observa- 
tionally determined LF bump is simply a result of a sys¬ 
tematic error introduced through the conversion from lu¬ 
minosity to absolute magnitude. The fact that the mag¬ 


nitude difference at FDU is not as pronounced as for 
the LF bump does not rule out a conversion problem. 
The bolometric corrections are functions of metallicity, Tefj 
and surface gravity and thus vary throughout evolution. 
ISalaris. Cassisi fc Weisd (l2002l ) have found that by changing 
the model atmosphere sets, bolometric corrections can differ 
by up to 0.1 magnitudes. This is not enough to account for 
the « 0.3 Mv offset between the models and observations 
at FDU. Even if we were to apply bolometric corrections 
from an empirically calibrated set of m odel atmospheres 
(e.g., iHoudashelt. Bell fc Sweigard 12000 ) . it does not rule 
out the presence of systematic errors in our conversion. To 
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reduce this source of error we us e independ e ntly c onverted 
visual magnitudes calculated by iLind et ^ (l2008l l in their 
study of NGC 6397. Their method, similar to ours described 
above, converts visual magnitude to luminosity by apply¬ 
ing a calibration from lAlonso, Arribas fc Martmez-Rogen 
lll999h . This calibration, too, is a function of (independently 
determined) metallicity and T^s. Magnitude was converted 
to luminosty through a 13.5 Gyr isochrone for the cluster 
l|R,ichard. Michaud fc Richerll2005l . which placed their stars 
in the mass range 0.78 — 0.79 M p,; similar to the 0.8 Mq 
that we have modelled here). The lhind et al.l (l2008l ) results 
are presented in the top right panel of Figure [S] Our stellar 
model is again denoted by the solid lavender curve and in 
this case, the luminosity is calculated directly from the equa¬ 
tions of stellar structure. The left (abundances as a function 
of Mv) and right (abundance as a function of luminosity) 
panels look remarkably similar with the expected factor of 
2.5 difference in the respective brightness scales. Mucciarelli 
(private communication) has also provided us with indepen¬ 
dently determined luminosities for M4 (Figure [3]). The two 
brightness scales {Mv and luminosity) yield the same be¬ 
haviour in each cluster, thus systematic uncertainties in our 
conversion are unlikely responsible for the magnitude dis¬ 
crepancy between theory and observation. 

We also have no reason to believe that there is a sys¬ 
tematic issue with bolometric corrections determined from 
model atmospheres at low metallicity. In fact one would 
expect the opposite to be true; bolometric corrections at 
high metallicity should be more uncertain. Synthetic colors 
perform quite well in the metal-poor regime, without sig¬ 
nificant deviations between metallicities of [Fe/H] = —1.5 
and [Fe/H] = —2.5 (Gasagrande, private communication). 
As the metallicity increases, however, the choice of atomic 
line-lists becomes increasingly important. We note that mi¬ 
croturbulence in low-metallicity model atmospheres is one 
possible source of systematic differences. Microturbulence 
can affect UV/blue wavelengths so if the velocities were to 
change in the low-metallicity regime, then the derived cor¬ 
rections would be systemically offset. As no evidence for 
such behaviour exists, a closer look at the physics of stellar 
modelling is required. 


4 ANALYSIS OF THE PHOTOMETRIC DATA 

The results from NGC 6397 and M4 raise a series of inter¬ 
esting questions that only a larger sample of GCs will help 
answer: 

(i) Do stellar models reproduce the structure of high- 
metallicity stars better than low-metallicity stars? 

(ii) Do uncertainties in the distance modulus make com¬ 
parison between theoretical and observed FDU and bump 
magnitudes (direct comparison hereinafter) too inconsistent 
to be meaningful? 

Unlike FDU which can only be identified by spectroscopic 
determination of abundances, the LF bump is readily iden¬ 
tified through the CMD and hence has been observed in 
many more clusters. By switching to photometric data a 
second evolutionary indicator in the main-sequence turn-off 
(MSTO) magnitude can also be employed as a check on the 
models. 


Comparisons between models and photometric GC data 
are common in the literature. There are three well tried 
methods by which we can compare theoretical predictions 
with empirical measurements of the LF bump: 


(i) By comparing the parameter = Vbump — Vhb; 

the V magnitude difference between the RGB bump 
and the horizontal branch at the RR L yr ae instability 
strip magnitude llFusi Pecci et al.l Il990l: IZoccali et al.l 


19991: ISalaris. Cassisi fc Weisd I 2 OO 2 I: iRiello et al.l | 20^ 


BjoiL^_Qiaba^ 20061 : 
Di Gecco et al J|2Q10|). 


iMeissner fc WeisX 


(ii) By comparing the parameter = 

I4isto — Vbump; the V mag nitude difference between 
the RGB bump and the MSTO dCassisi et al.llioilh . 


(iii) By comparing the photometrically determined abso¬ 
lute mag nitude of the LF bump to that predicted by the 
mode ls llKing.JD<i_Cpsta^JM)CTnarqu ^ [l9^; j Alongi et al.l 


- - )nm 

I 1991 I : iBiork fc Chaboverl I 2 OO 6 I : lAngelou et al.l 12012!) . Each 
stellar model requires a bolometric correction to convert lu¬ 
minosity to Mv whilst the distance modulus is required to 
convert the observed visual magnitude to absolute magni¬ 
tude. 

The hrst two methods have the advantage of being in- 
dependent o f dista nce and reddening. As discussed by 
ICassisi et al.l (l201ll L the first method does introduce some 
uncertainty with respect to the placement of the observed 
HB level for GGs with blue HB morphologies and in theo¬ 
retical predictions of the HB luminosity (dependent on each 
code and their predicted ‘‘He core mass at the ^He ignition at 
the RGB tip). The third method introduces an extra source 
of uncertainty because it relies on accurate determination of 
the distance modulus (as well as suffering from uncertainty 
due to interstellar reddening). In order to understand the 
level of error introduced from uncertainties in the distance 
modulus we employ both the second and the third method 
in our investigation of the RGB stellar models. 

The natural method for comparing the key theoreti¬ 
cal and photometric indicators of stellar evolution has been 
through fitting of isochrones to the CMD. In this study we 
calculate individual stellar models at distinct metallicities 
and fit third order polynomials to create what are essen¬ 
tially coarse isochrones. Table U outlines our sampling of 
the GC metallicity range and indicates the corresponding 
initial masses for each metallicity that yield a MSTO age of 
12 Gyr. 

We note that it is common to require isochrones with 
ages greater than the age of the Universe to match the 
luminosity of the L F bump (see iRiello e t al.l l2003l and 
ICassisi et al.ll2ojTl l . As ICassisi et al.l ( 2 OIIII comprehensively 
discuss, the reason that such old ages are needed is be¬ 
cause the underlying stellar models fail to reproduce the 
RGB bump brightness for an age appropriate to the GC as 
measured from its MSTO brightness. Thus there exists clear 
motivation to focus on the underlying stellar models and the 
included physics. 

GCs are observed to have enhanced a-element abun¬ 
dances compared to scaled-solar values and their metallic¬ 
ities are often expressed in terms of the total metallicity , 
[M/Hj. This is defined bv ISalaris. Chieffi fc Stranierd (1199311 
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z 

[M/H] 

Mass (Mq) 

0.02 

0.17 

1.017 

0.011 

-0.09 

0.956 

0.006 

-0.35 

0.900 

0.0036 

-0.57 

0.854 

0.002 

-0.82 

0.828 

0.0011 

-1.09 

0.812 

0.0006 

-1.35 

0.805 

0.00035 

-1.58 

0.800 

0.0002 

-1.83 

0.800 

0.00011 

-2.09 

0.800 

0.000063 

-2.33 

0.800 


Table 4. The metallicities of the grid of models used to generate 
isochrones in Figures [4| and [5] The third column indicates the 
initial mass of the star representative of the stellar structures 
populating the RGB at an age of 12Gyr. 


as 

[M/H] = [Fe/H] + log(0.638 x + 0.362). (5) 

This metallicity definition has been used in similar studies 
and is adopted here. 

As we have mentioned, this type of comparison has 
been carried out extensively in the literature. In t heir analy¬ 
sis of the AFg™*’ parameter, iRiello et al.l ll2003ll find good 
agreement between theory and observations at higher metal¬ 
licities but note significant discrepancies at low metallic- 
ity. Investigations of l ow- metallicity GCs include those by 
iFusi Peed h al.l lll990l ~) and iDi Cecco et al. I ll20l^ , who also 
used the parameter to determine that at [M/H] 

< —1.7, models underesti mate the magn i tude of the LF 
bump by at least 0.4 mag. iDi Cecco et al.l (l2010f l also con¬ 
ducted tests on the effects of the microphysics and deter¬ 
mined that models with a-element and CNO-enhancements 
could not account for the discre pancy nor could th e re¬ 
vised solar heavy-element mixture. ICassisi et al.l (1201 il l, who 
employed the parameter, found that the theoret¬ 

ical bump was too bright by 0.2 magnitudes on average, 
but discrepancies of AMy ^ 0.4 magnitudes in the lowest- 
metallicity clusters were within the uncertainties. In their 
test of the stellar microphysics, iBiork fc Chabo\^ (l2006l f 
have determined that the uncertainty in theoretical values 
for the LF bump magnitude varies with metallicity between 
+0.13 and —0.12 magnitudes at [Fe/H] = —2.4 and between 
+0.23 and —0.21 magnitudes at [Fe/H] = —1.0. The dom¬ 
inant sources of uncertainty were attributed to a-element 
abundance, the mixing length parameter, and the high- 
temperature opacities, all of which are increasingly impor¬ 
tant at higher metallicity. We stress their main result is that 
from a purely theoretical perspective, there is more scope 
for unce rtainty in the meta l-rich models. We also note the 
work bv lTroisi et al.l (l201ll ) who used the empirical bright¬ 
ness difference between the LF bump and the point on the 
main sequence that is at the same colour as the bump. Their 
models also underestimated the magnitude of the LF bump 
unless they employed an initial He mass fraction Y = 0.2 
which is lower than the Big Bang Nucleosynthesis value. 


4.1 Comparison Across the Globular Cluster 
Metallicity Distribution 

4-1.1 Results and Discussion 

In Figure [4] we plot the magnitude of two key evolutionary 
features as a function of metallicity. Matching the MSTO 
luminosity (panel 2^) and LF bump luminosity (panel UJo) 
serve as initial tests for our models. To remove the effects of 
distance and reddening uncertainties, we also compare the 
magnitude difference between the locations of the MSTO 
and LF bump (panel |4};). In these three panels the observa- 
tional data are pl otted with black circles and taken from 
iNataf et al.l (l2013fl who performed statistical analyses on 
HST photometry. Their work yielded the MSTO magnitudes 
for 55 clusters to which we compare our models in panels 2^ 
and 2}:- These 55 clusters are a subset of the 72 systems 
shown in 2t for which they determined the LF bump mag¬ 
nitudes. 

We have constructed a grid of models that vary in 
mass and metallicity to coincide with MSTO ages of 8 Gyr 
(lavender curve), 12 Gyr (blue curve) and 15 Gyr (vermilion 
curve). We fit third order polynomials through the quanti¬ 
ties of interest to interpolate across the entire globular clus¬ 
ter metallicity range. The metallicity spacing of our grid is 
listed in Table 21 and is the same for each isochrone. The 12 
Gyr isochrone will serve as somewhat of a fiducial track in 
our analysis and in Table 21 we indicate the corresponding 
mass at each metallicity for this case. In all models we as¬ 
sume an initial hydrogen abundance of 0.75 and we assume 
a solar scaled ASG09 metallicity mixture with no a-element 
enhancement for the initial composition and opacity mix¬ 
tures. 

To quantify the discrepancy between theory and obser¬ 
vation at the LF bump, we compare our models with the 
lines of best fit to the data in pan els Hb and2t (f^., equa¬ 
tions 1 and 2 in iNataf et al.ll^013l f 

Mv.bump = 0.600 + 0.737([M/H] + 1.110) (6) 

AMv = 3.565 - 0.549([M/H] + 1.152). (7) 

These fits are translated to give upper and lower envelopes 
to the data in order to provide an indicative uncertainty 
at each metallicity. The discrepancy between each of these 
three fits and the 12 Gyr isochrone can be found in Panel 
2Ji, (top) for the direct comparison method and in panel 21^ 
(bottom) for the distance independent parameter. 

In panel 2K we run into the familiar result (see also 
figure 2 of Cassisi et al. 2011) that 15 Gyr isochrones are re¬ 
quired to match the MSTO in many GGs. It is expected that 
given the age constraints of Bi g Bang Cosmology and the 
cosm ic microwave background (IPlanck Gollaboration et al.l 
l2014h that the 12 Gyr isochrone should provide an upper age 
limit to the clusters. The data demonstrates a clear spread 
in age at each metalliciy with a collection of clusters cen¬ 
tred around [M/H] « —0.3 dex well modelled by the 8 Gyr 
isochrone. 

The stellar models also tend to predict LF bump mag¬ 
nitudes that are too bright compared to the observational 
data. There is some agreement at the metal rich end of the 
spectrum above a metallicity of [M/H]= —0.8. At the low¬ 
est metallicites the models differ significantly - by up to 0.8 
magnitudes for some clusters (see panel 2Ji)- The fact that 
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Panel a Panel b 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 


[M/H] 


[M/H] 


Figure 4. Panel a: MSTO magnitudes for 55 GCs (black circles with associated uncertainties. iNataf et ^1201.111 as a function of [M/H]. 
We include predictions from our stellar models assuming a MSTO age of 8 Gyr (lavender), 12 Gyr (blue) and 15 Gyr (vermilion). Panel 
b: The empirically determined location of the LF bump in 72 GGs as a function of [M/H]. The data are from th e sources de s cribed 
in panel a and retain their symbols. The grey curve in this panel is the line of best fit to the data determined bv iNataf et al.l ll2013l l 
and described by Equation 6. Panel c: The distance independent parameter AMy as a function of [M/H]. The data are from 

the sources descri bed in panel a and retain their symbols. The grey curve in this panel is the line of best fit to the data determined by 
iNataf et al.l ll2013f ) and described by Equation 7. Panel d; The Euclidean distance between the LF bump magnitude determined from the 
12 Gyr isochrone and respective best, lower bound and upper bound fits to the observational data in panel b. Panel e: The Euclidean 
distance between the AMy parameter determined from the 12 Gyr isochrone and respective best, lower bound and upper bound 

fits to the observational data in panel c. 


our models fail to reproduce the slope of the observational 
data suggests that we need to give careful consideration to 
the internal stellar processes and how they operate across 
different metallicities. 

Similarly, we find some agreement for the distance- 
independent AMv^^jP parameter at the highest metal¬ 
licities. The level of discrepancy outlined in panel |4^ is 
consistent with a roughly constant offset of SMy « 0.4. 


This is s lightly highe r than the value of SMy « 0.2 
found by ICassisi et al.l (l201ll ) but their models better re¬ 
flect GC abundances (see below for discussions on "^He and 
a-elements). It is essential to point out the ‘better’ agree¬ 
ment implied by the 8 Gyr isochrone. This track significantly 
overestimates the MSTO luminosity and predicts an incor¬ 
rect age in the low-metallicity clusters but in doing so, better 
matches the AMy value. 
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-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 

[M/H] 


Figure 5. Same as Figures |4^ and |4l> except assuming 3% in¬ 
crease to the distance modulus for each cluster. Note the addi¬ 
tional results for age= 4.0 Gyr and Y= 0.40. 


4-1.2 The Distance Modulus and Metallicity Scale 

Our results, irrespective of how the comparison is carried 
out, imply that the stellar models do not correctly repro¬ 
duce the stellar structure at low metallicity. The MSTO and 
LF bump magnitudes are underestimated (too bright) in our 
models and we do not match the magnitude difference be¬ 
tween these two evolutionary indicators. Naively, if we are 
unable to match the earlier evolutionary events then we have 
no right to expect to match the later events, but it is pos¬ 
sible to shift the magnitude of the LF bump (essentially) 
independently of the MSTO. The LF bump is of course de¬ 
pendent on the depth of FDU: if FDU is deeper, then the 
hydrogen-burning shell will encounter the composition dis¬ 
continuity at a lower luminosity. Improvements to stellar 
models can therefore be attained through either refinement 
of the microphysics or changes in the treatment of mixing. 

Before we investigate the role of the physics included 
in the stellar models, we first determine whether uncertain¬ 
ties in the observational data can help reduce the disagree¬ 
ment with theoretical predictions. One of our key ques¬ 
tions pertains to the usefulness of direct comparisons and 
the uncertainty surrounding the distance modulus. A range 
of methods are used to determin e the d istanc e to GCs. 
The distance moduli listed in the iHarrisI lll996ll catalogue 
(2010 edition) are predominately determined from c alibrat¬ 
ing th e mean V magnitude of the horizontal branch. iHarrisI 
(Il996l l cites an uncertainty in determining this mean mag¬ 
nitude of at least ±0.1 Mv- There are also unquantified dif¬ 


ficulties i n the magnitude-meta llicity relation used for cal¬ 
ibration (iDemaraue et al]|200Clh . In the case of NGC 6397, 
the blue horizontal branch introduces much uncertainty for 
this method. Its distanc e modulus was deri ved by fitting 
isochrones to the MSTO (lAlcaino et al.|[T987^ . In fact most 
of these techniques used to determine the distance modulus 
require some form of calibration from isochrones. 

Figure [5] demonstrates how a systematic offset in the 
distance modulus determination would help to reconcile 
models with the observational data. We have reproduced 
panels UK and UK but multiplied the distance modulus for 
each cluster by an ad hoc factor of 1.03. The 12 Gyr MSTO 
models now matches the upper extent of the MSTO mag¬ 
nitudes across all metallicities. Directly matching the LF 
bump with our isochrones requires that some metal-rich 
clusters possess a '^He content of up to Y=0.4, however, this 
isochrone is unable to simultaneously provide a lower limit 
to the MSTO magnitude in these systems. Further changes 
to the input physics are seemingly necessary to account for 
the systems with the brightest LF bumps. In fact, 4 Gyr is 
probably not a realistic age estimate for a GG given their 
metallicity. Nevertheless, it is somewhat pleasing that pa¬ 
rameters within known observational limits can reproduce 
most of the data; a more reasonable state of affairs than re¬ 
quiring stars with ages greater than the Hubble Time. The 
uncertainties introduced by the absence of boutique mod¬ 
elling (where all known cluster parameters are considered) 
and self-consistent determination of the distance modulus 
(where fitting is done with the code generating the models) 
may contribute to the need for this factor 1.03. 

The chosen metallicity scale is also a source of un- 
certainty and has been discussed in previous studie s 


dRiello et al.ll200)j^ Di Cecco et al.llioi^ : ICassisi et al.ll20li] ~) . 


iDi Gecco et al.l l 20ld L in particular, find that by adopting 
two different sets of independently determined metallicites 
for their GG sample, their conclusions remain unchanged. 
The models still underestimate the bump magnitude. We 
note that a systematic offset of about 0.5 dex in [M/H] would 
help reduce some of the inconsistency between theory and 
observation but there is currently no evidence for advocat¬ 
ing this. Ghanges to the distance modulus and metallicity 
scale are equivalent to applying luminosity translations in 
the HR diagram. So whilst ‘corrections’ to the observational 
data may assist with the direct comparison methods they do 
not help with the distance independent comparisons. The 
AMv bu^i? parameter is set by the stellar physics and es¬ 
sentially a constant offset at low metallicity. 


5 INVESTIGATION OF THE STELLAR 
PHYSICS 

The choice of stellar physics has some role to play in rec¬ 
onciling the differences between theoretical and observa¬ 
tional determinations of the MSTO and LF bump magni¬ 
tudes and, in particular, reducing the luminosity difference 
between these events. In the following sections we quantify 
how choices in the modelling affect the bump magnitude and 
the AMv bimii? parameter. Where possible, we generate new 
isochrones to illustrate the effect of the stellar physics on our 
parameters of interest, but in some cases we compare indi¬ 
vidual stellar tracks because it makes the analysis simpler. 
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Figure 6. Top panel: stellar tracks for stars representative of 
GC PIE populations (see definition in text) at a metallicity of 
Z=0.00011. Mass and helium abundances are as specified in the 
legend. Each model reaches the MSTO with an age of 12 Gyr. 
Bottom panel: penetration of the convective envelope as a func¬ 
tion of luminosity for the respective models. Dotted vertical lines 
mark the bump location in each case. 


We note that because of the short lifetimes of the SGB and 
RGB, single evolution tracks coincide almost exactly with 
isochrones. Our models that turn off the main sequence at 
an age of 12 Gyr will serve as an appropriate proxy for the 
isochrone and once the effects of the microphysics are quan¬ 
tified, they can later be incorporated in the grids of models 
that generate the isochrones. 


and pristine cluster gas determines whether the second gen¬ 
erations are ‘extreme’ population stars (sometimes formed 
from pure ejecta) or belong to the ‘intermediate’ popula¬ 
tion that give rise to the o bserved anticorrelati ons (mixed 
ejecta and cluster gas, see ICarretta et al.l l2009l for discus¬ 
sions on [PJrimordial, [I]ntermidiate and [E]xtreme popula¬ 
tions [PIE]). The extreme population is not evident in ev¬ 
ery GG, however, it is predicted that all second generation 
stars are necessarily enhanced in helium compared to their 
primordial counterparts. A spread in '^He abundance also 
straightforwardly ex plains the horizontal branch morphol¬ 
ogy i n many clusters IP’Antona fc Cai^l2008l : iMarino et al.l 
I 2 OI 4 I and is consis t ent with observed split main sequences 
i Piotto et al.ll2007l : iMilone et al.l[201^ . 

The effect of ^He on the stellar models is well 
documented in the litera t ure - originally in a series 
of papers by Iben Jlben|_[l968|;_llben_&_^kulkneil Il96^ 
and more recently b y lGassjsi ^_Salarjsfc Pietoinfernil (l2013l) . 
ISalaris et al.l ll2006l ) and iGassisi et al~ (2011) place an up¬ 
per limit of AA('^He) = 0.05 on the width of the main 
sequence for those clusters that do not display split evo¬ 
lutionary tracks. This in turn decreases the TO age by no 
more than « 0.5 Gyr. In clusters such as u Gen, a spread in 
metallicity and split main-sequences imply a wider span in 
age, however, these complex systems are in the minority. In 
most cases, GGs are dominated by second generation stars 
with an inferred separation in age of at most a few hundred 
Myr (iGarretta et al.l2010bl) . Stellar models representative of 
these stars only further exacerbate the discrepancies demon¬ 
strated in Figure m In Figure [S] we plot the HR-diagram (top 
panel) and penetration of the convective envelope (bottom 
panel) for three stellar models representative of the PIE pop¬ 
ulations at a metallicity of Z=0.00011. Each model, with its 
combination of initial mass and '*He content leaves the main 
sequence with an age close to 12 Gyr. Our models are by 
design the same age and metallicity so there are minor dif¬ 
ferences in their MSTO luminosity but they differ in tem¬ 
perature. In the bottom panel the luminosity of the bump 
for each model is marked by the dotted vertical lines. The 
impact of enhanced helium on the envelope opacity reduces 
the extent of FDU, and delays the onset of the bump (« 
0.2 Log L/Lq). The second generat ion models (see also the 
Y=0.40 isochrone in Figure and [Claris et ^ I 2 OO 6 I ) are 
brighter, contrary to the need for fainter MSTO and LF 
bump magnitudes demons trated in Figure 2] It follows that, 
as per iTroisi et al.l (l201ll) . primordial '^He abundances be¬ 
low the accepted Big Bang nucleosynthesis level would allow 
models to better reproduce observations. 


5.1 Mass and ^He-Enhancement 

The existen ce of multiple populations in GGs is well docu¬ 
mented (see lGratton. Garretta fc Bragagliall2012l and refer¬ 
ences therein). The constituent stars divide into chemically 
distinct groups that are best explained by explicit stellar 
generations. The most commonly proposed scenario is that 
a fraction of first generation stars pollute the environment 
from which the subsequent generation(s) form (which we 
herein refer to as second generation). More-massive stars in 
the first generation, having experienced the advanced phases 
of stellar nucleosynthesis, burn hydrogen via the CNO cy¬ 
cles. The degree of mixing between first-generation ejecta 


5.2 Opacity 

The models presented thus far do not account for the ef¬ 
fect of a-element enhancement on the evolution/isochrones. 
The inclusion of a-elements has a two-fold effect. First, it 
increases the amount of available GNO material which im¬ 
pacts upon the nuclear burning. Second, with a larger reser¬ 
voir of metals, there is a contribution to the opacity that 
is predominantly manifested as a reduction of the effective 
temperature. As a function of [M/H], the inclusion of a- 
element enhancement results in only minor changes to the 
MSTO luminosity. 

In the top panel of Figure [3 we plot the bump mag- 
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Figure 7. Top panel: The predicted bump magnitude for 12 Gyr 
isochrones with various levels of a-element enhancement. Bottom 
panel: HR diagram for the 12 Gyr isochrone with AGS09 solar 
scaled oxygen abundance but different selections of a-element en¬ 
hancement in the heavy mixture opacity tables. 


nitude as a function of [M/H] determined from 12 Gyr 
isochrones with [a/Fe]= 0.0 (blue curve), [a/Fe]= 0.2 (ver¬ 
milion curve) and [Q/Fe]= 0.4 (lavender curve). The heavy 
element mixtures for the respective opacity tables can be 
found in the appendix, a-element enhancement shifts the 
bump to brighter magnitudes in the more metal-rich mod¬ 
els. This effect is on account of the increased O available 
in the nuclear burning chains. In the scaled-solar case a- 
element enhancement of [Q/Fe]= 0.4 shifts the bump by 
approximately 0.15 mag. 

The bottom panel of [7] illustrates the effect of the dif¬ 
ferent opacity tables on the individual stellar models. We 
have calculated evolutionary tracks based on a star with 
metallicity Z=0.02 and TO age of 12 Gyr used in generating 
our isochrones. In these calculations we employ a-enhanced 
opacity tables but do not change the composition of the 
models. The CNO abundances in our energy generation net¬ 
work are kept at the ASG09 scaled-solar values, which re¬ 
moves the impact of increased CNO material on the bump 
location. The opacity tables used are described in the top 
panel and detailed in the appendix. We also include a calcu¬ 
lation (black curve) with the GN93 scaled-solar composition 
and find variations of up to 0.04 Log L/Lq are possible in 
the bump location. 


5.3 Mixing Length Parameter 


The free parameter in the MLT is also a source of un¬ 
certainty in the stellar models. q.mlt is calibrated to 
the Sun, a star halfway through its main-sequence life¬ 
time which possesses a shallow convective envelope (in 
mass). This value is then applied across the stellar mass 
and metallicity dist ribution and through a ll pha ses of 
evolution. Work by iMagic. Weiss fc Asplundl ll2015lf sug¬ 
gests that this should in fact vary with evolutionary 
status but the solar calibrated mixing length parameter 
does allow a majority of the mo del libraries e.g., BA STI 
(|Pietrinferiij_et_^ 2004ll. Padova fB ressan et al.l]2( )12lf and 


Victoria I VandenBerg. Bergbusch fc Dowl^ 200611 to match 

the effe ctive temperature of RG B stars in GCs (see also the 
work bv ISalaris fc Cassisill20ld l. The value of umlt adopted 
will significantly shift the model’s position in the HR dia¬ 
gram because the convective efficiency will change the sur¬ 
face temperature. The effect is greatest on the RGB where 
the stars possess deep convective envelopes. Although not il¬ 
lustrated here, we find that the depth of FPU is not altered 
significantly by this parameter (see also lAlongi et al1ll99l[l 
and we find variation of approximately 0.05 mag between 
ctMLT = 1.60 and ocmlt = 1.85. 


5.4 Equation of State 


MONSTAR typica l ly em ploys the fitting formula by 
iBeaudet fc Tassoull (Il97llf whilst partially-ionized regions 
are treated wit h the Saha equation as described by 
iBaerentzenI (Il965fl . In the event of convergence issues, a com¬ 
putationally more expensive analytic solution can be calcu¬ 
lated but this method still requires the numerical evaluation 
of the Fermi-Dirac integrals. To aid in this investigation we 
have a dded the OPAL EOS tables from [Rogers fc NavfonovI 
ll2002ll . We also include models from MONSTAR in¬ 
corporating the Timmes EOS (iTimmes &: Arnett 1999h 
and the Helmholtz equation of state ll Timmes fc Swesra 
I 2 OOOII which is the tabulated form of the Timmes EOS 
i Gonstantino et al The Timmes EOS has allowances 

for “simple ionisation” that includes a simple two-level 
hydrogen-like atom model for ionization. In cases where 
equations of state are blended, a linear transition occurs 
over the temperature range T = 2 — 2.5 MK. 

We have calculated the evolution of a star with mass 
M= 0.8 Mq and metallicity Z=0.00011 using five differ¬ 
ent EOS combinations (Figure|8]). Employing the Helmholtz 
EOS in the high temperature regime produces a systemati¬ 
cally higher at a given luminosity compared with other 
choices for the EOS. The depth of FP U for each combina¬ 
tion is plotted in the adjacent panel. The IBeaudet fc Tassoull 
1 I 1971 II fitting formula, which is the default choice in MON¬ 
STAR, provides the shallowest penetration of the convective 
envelope of all five EOS combinations. This choice of EOS 
therefore yields the greatest discrepancy with the empiri¬ 
cally determined LF bump magnitude. The deepest pene¬ 
tration of the convective envelope is achieved by using the 
OPAL EOS and leads to approximately 0.1 magnitude bet¬ 
ter agreement with the empirical data. 
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Figure 8. Left panel: Evolution through the HR-diagram of the stellar model with metallicity Z=0.00011. Each curve corresponds to a 
different choice for the EOS. Right panel: Penetration of the convective envelope for each EOS choice. 




[M/H] [M/H] 


Figure 9. Left panel: LF bump magnitudes as a function of metallicity from the observations specified in Figure |4] (black circles). 
The coloured curves denote our 12 Gyr isochrone with various levels of overshoot at the base of the convective envelope. Right panel: 

parameter for the observational data specified in Figure [4] (black circles) and for our 12 Gyr isochrone with various levels of 
overshoot (coloured curves) at the base of the convective envelope. 


5.5 Overshoot 

Once it had been identi fied that stellar models overestimate 
the bnmp Inminosity | King. Da Costa fc Demarauel I198 !tI : 
iFusi Peed et al.lll990h . Alongi et al.l (ll99ll ) suggested that 
overshooting by 0.7 — 1.0 pressure scale heights (77p) at the 
base of the convective envelope would provide the required 
shift in magnitude to reconcile theory with observations. 
Using the formalism outlined in iJ5] we find that values of 
fos = 0.05 — 0.1 significantly improve the agreement with 
the observations - both with the direct comparison of the 
bump magnitude and with the Mv parameter. In Fig¬ 

ure [9] we plot the 12 Gyr isochrone with no overshoot (blue 
curve), fos = 0.05 (lavender curve) and fos = 0.1 (vermilion 


curve). We do not expect that all GCs are 12 Gyr old, but if 
these tracks can provide appropriate extrema to the data in 
Figure [9] then a spread in age and overshoot efficiency may 
allow models to reproduce the observations. 

The isochrones reproduce the gradient of the Mv bimu? 
parameter, but do not match the bump magnitudes. The 
isochrones still do not provide an upper limit to the some 
of the most metal-poor clusters but do satisfactorily explain 
intermediate-metallicity and metal-rich clusters. This may 
hint that at the need for an overshoot formalism that is a 
function of metallicity such that: 

• for [M/H] > 0 no overshoot is required; 
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• between —0.8 < [M/H] < 0.0 overshoot increases to 
fos « 0.05; 

• between —1.3 < [M/H] < —0.8 overshoot increases to 

fos « 0.1; 

• for [M/H] < —1.3 overshoot fos > 0.1 is required. 

It is interesting to note that our models with overshoot 
cannot simultaneously provide a lower limit to the absolute 
magnitude of the bump and the My parameter. In the 

left panel of Figure |9] the most efficient form of overshoot 
we have modelled predicts a bump location that is approx¬ 
imately 0.3 magnitudes too bright in the most metal-poor 
clusters. If we include more efficient overshoot such that our 
curve encompasses the bump magnitude of these systems, 
it will necessary reduce the separation between the MSTO 
and LF bump. This is in spite of the fact the vermilion 
curve {fos = 0.1) currently forms a satisfactory lower enve¬ 
lope. The discrepancy between the two panels implies that 
the models predict a MSTO magnitude that is too bright for 
a 12 Gyr isochrone - an outcome of different input physics. 
Such a result reflects the fact 15 Gyr isochrones are required 
to account for the MSTO luminosities in Figure [4^. 


5.6 Other Input Physics 


Determinations of the ke y hydrogen burning reaction rates 
have changed little since [Harris et al.l (|l983h save for the 
significant revision of the ^^N(p, 7)^®0 rate which we have 
discussed in and use in the current calculations. One 
possible source of improvement lies with the recent mea¬ 
surement of th e iron opacity at te mperatures pertaining to 
solar interiors. iBailev et al.l (l2015h find that the measured 
iron opacity is much higher (30- 400%) than previously pre¬ 
dicted . In their series of tests ISalaris. Chieffi fc Stranierd 
lll993tl demonstrate the significant role that the high tem¬ 
perature opacities have on the luminosity of the MSTO. An 
increase in Fe opacity will help reduce the MSTO luminos¬ 
ity. Only a small improvement is likely because the change is 
wavelength dependent and Fe is just one of many elements 
that contribute to the opacity. 

The inclusion of additional physical processes is an 
obvious avenue of pursuit with the role of atomic dif¬ 
fusio n especially a contentious issue (see iGruvters et al.l 
I 2 OI 3 I and references therein). If the process is not inhib¬ 
ited then atomic diffusion will re duce the MS lifetim e 
and hence the turn-off luminosit y (IChabover et al.ll200iri . 
Observational evidence however i Gratton et aJ 200ll~) sug- 
lests that the process is indeed inhibited. ICassisi et al.l 
2 OIII I state that the overall effect from including diffusion 
on Mv bimfi? is a reduction of approximately 0.05 mag. 
The role of rotation and magnetic fields are also uncer¬ 
tain but their inclusion will have consequences for t he en - 
tire evolution. Chaname, Pi nsonneault fc Terndrud (l2005l 'l 


and jPalaci^ et al.l ( 200d l demonstrate how the treatment 


of angular momentum impacts upon the stellar models. 
In their comparison of different transport prescriptions, 
[Palacios et al.l ([2006[ ) found a negligible difference in MSTO 
luminosities and the depth of FDU. However, each of the 
rotating models still predicted a higher LF bump lumi¬ 
nosity than in the non-rotating models. The higher-mass 
higher-metallicity (M > 1.25 M©, Z models by 

[Charbonnel fc Lagardel ( 201Clh . on the other hand, demon¬ 



Figure 10. Top panel: A(Li) as a function of magnitude for NGC 
6397 ([Fe/H] = —2.1). We include appropriate models with three 
variations of the thermohaline mixing parameter Ct ■ Blue curves 
denote Ct = 120, lavender denote Ct = 300 and vermilion curves 
denote Ct = 1000. Bottom Panel: Panel b: [C/Fe] as a function 
of luminosity in NGC 6397 (black circles) and NGC 5466 ([Fe/H] 
= —2.2, grey circles). Model data are as described for the above 
panel. Sources for the observational data (in both panels) are 
described in the main text. 


strate a significantly lower bump luminosity in their rotating 
models compared to their non-rotating models. 


6 SECOND ANALYSIS OF THE 

SPECTROSCOPIC DATA: THERMOHALINE 
MIXING 

In this section, we compare observations and theoretical 
predictions of the evolution of two species that are in¬ 
volved in different channels of hydrogen burning: the fragile 
lithium and the more robust carbon. We have run a grid of 
models rather than conducting boutique modelling of each 
cluster. All models were run with initial X(H)= 0.75 and 
X(Li)= 9.39 X with an assumed TO age of 12 Gyr. 

Overshoot (with fos = 0.075) has been included and ensures 
that the parameter is consistent with the value 

determined from the respective CMDs (but does not neces¬ 
sarily reproduce the bump magnitude). In order to test the 
efficacy of the thermohaline process, we have applied artifi¬ 
cial offsets in magnitude to the models so that the beginning 
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Figure 11. A(Li)as a function of magnitude for six of the clusters described in Table [2] For each panel we calculate stellar models 
that best represent the metallicity of the respective clusters. We include three variations of the thermohaline mixing parameter Ct as 
described in the legend. We also indicate the location of the bump magnitude with the vertical dashed line and the sources from which 
the observational data are taken. The error bar in the top-left panel indicates the typical uncertainty in the data. 


of extra mixing coincides with the photometrically derived 
LF bump. 

In the top panel of Figure[T0]we once again plot lithium 
observations of NGC 6397 (Lind et al. 2009, black circles) 
but include models with three different choices for the free 
parameter in the diffusive-thermohaline mixing theory (Ct). 
Although Ct ~ 1000 has been the preferred value in the lit¬ 
erature, matching abundances in many stellar environments, 
values between Ct = 120 — 300 give a far better match to 
Li in this system. The same models are included in the bot¬ 
tom panel of Figure fTOl in our comparison to [C/Fe] data in 
NGC 6397 and NGC 546 6. [C/ Fe] data for NGC 6397 were 
taken from iBrilev et al.l (Il990ll but determined by assum¬ 
ing that there is no oxygen over- abundance in that sys tem. 
[C/Fe] in NGC 5466 is taken from IShetrone et al.l (l2010l ). We 
expect that the models are representative of both clusters 
despite the small difference in metallicity. NGC 5466 was 
one of the three metal-poor GCs studied bv lAngelou et all 


(l2012li in their investigation of thermohaline mixing. Of the 
three, this was the only system in which thermohaline mix¬ 
ing appeared to reproduce the depletion of [C/Fe] along the 
RGB. NGC 6397 seems to display the same ambiguity as 
its more massive counterparts: it is difficult to identify the 
onset of extra mixing from the [C/Fe] data and one could 
argue that it begins well before the expected bump. [C/Fe] 
in NGC 5466, having been homogeneously analysed from a 
cluster of comparable size to NGC 6397, serves as the best 
analogue to the Li data. 

We find that a mixing efficiency of Ct = 1000 depletes 
Li much too fast (top panel Figure fTOl) but is able to account 
for the depletion of [C/Fe] in both clusters (bottom panel 
Figure [Toll . In Figure ITT] we investigate lithium depletion the 
clusters listed in Table [2] with the exception of NGC 2808 
where the data does not include stars experiencing extra 
mixing. In most cases our models provide a lower limit to 
the post-FDU Li abundances (on account that we have not 
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tuned the initial abundances for each cluster). To guide the 
eye, we have artificially translated the Ct = 1000 models to 
a higher post-FDU abundance (dotted vermilion curves) to 
indicate how depletion would proceed with boutique abun¬ 
dance choices of each system. In nearly all GCs the models 
with Ct = 120 and Ct = 300 provide preferable fit to the 
Li data. The result is not entirely clear cut, especially in 
in the case of NGC 1904. However, we note that it is only 
the stars that have severely depleted Li that introduce some 
doubt. Those stars that are experiencing extra mixing (just 
after the bump) are better reproduced by the lower Ct val¬ 
ues thus these are our preferred values for reproducing the 
depletion of Li in GGs. In the bottom panel of Figure [TOl it is 
clear that with values of Ct = 120 and Ct = 300 not enough 
carbon is processed to match the trend of [C/Fe]. We are 
unable to simultaneously account for [C/Fe] and A(Li) with 
the same free parameter. Matching both these abundances 
poses a challenge for any postulated mixing mechanism. The 
fragility of Li and robustness of C must be reproduced. 


It may be that thermohaline mixing is not 
responsible for the surface composition changes 
of RGB stars. Objections to the mechanism are 
based on theoretical gr ounds and 3D hydrody- 
namical simul ations _ l[Tn3xjer j_Gnjaud fc Stellmrmd 


20 111 : 


Denissenkov fc Merrvfield 


20 111 : 


Brown. Garaud fc Stellmachl l2013ll that suggest the 


mixing is inefficient in the RGB regime. But we have de¬ 
scribed how the mechanism at work here can be considered 
a physically based phenomenological model. Rather than 
mixing to a constant temperature or mass location from 
the hydrogen burning shell, material is transported to 
where ^He burning creates an inversion in the /r profile. 
IChurch et all (120141 ) show how this location depends on the 


shell-burning conditions that change with RGB luminosity. 
In this extra mixing formalism we have varied the free 
parameter that controls the mixi ng rate. We note that 
IPenissenkov fc Pinsonneauld (l2008l l have advocated that a 
deeper (than the /j. inversion) but slower mixing may be 
required. It is not clear yet if this scenario could explain 
the abundance trends here. 


We have modelled thermohaline mixing as a diffusive 
process. Our evolution code employs a one-dimensional dif¬ 
fusion algorithm to model a process that is inherently three 
dimensional and advective. Gonvection is characterised by 
streams (or plumes) of upward and downward travelling ma¬ 
terial. The majority of the nuclear processing in the extra 
mixing region occurs where material turns over near the 
hydrogen-burning shell. A diffusive treatment of mixing as¬ 
sumes that the composition difference between convective 
elements at the same mass coordinate is negligible. If the 
mixing speed in the extra mixing region is much slower 
than convection (as it generally is believed to be), then the 
turnover time is closer to the reaction timescale. We there¬ 
fore expect a difference between the compositions of material 
approaching and returning from the envelope at the same 
mass coordinate. It may in fact be necessary to develop a 
more realistic treatment of convection to better model the 
extra mixing process. Godes that calculate mixing via ad¬ 
vective streams may prove to b e enlightening (ICannonll 19931 : 
IStoklllioO^ : iGhurch et al.ll200^ . 


7 CONCLUSIONS 

We have utilised photometric and spectroscopic measure¬ 
ments of GCs to investigate the current state of red giant 
branch modelling. Observations at key evolutionary epochs, 
namely the MSTO, FDU and LF bump, provide a series of 
constraints for our stellar models. Li determinations in seven 
GCs covering a factor of ten in metallicity have served as the 
key spectroscopic diagnostic. The surveys of NGC 6397, the 
most metal-poor cluster in our sample, and M4, the most 
metal-rich, are of particular interest as they cover a lumi¬ 
nosity range that extends from the subgiant branch to be¬ 
yond the LF bump. Homogeneous Li data, spanning the key 
evolutionary indicators has allowed us to confirm the long- 
suspected ambiguity associated with inferring the onset of 
RGB mixing from [C/Fe]. The many uncertainties in de¬ 
termining the [C/Fe] combine to give the appearance that 
mixing begins prior to the LF bump. In the seven clusters 
analysed here the onset of lithium depletion coincides with 
the photometric bump magnitude. 

We determined that the abundance changes predicted 
by FDU in the models are consistent with spectroscopic 
measurements, but we also found evidence that the luminos¬ 
ity at which mixing begins in the low-metallicity clusters is 
slightly overestimated. Data from a larger sample of clusters 
and predictions from other stellar evolution codes will help 
confirm this behaviour. The luminosity that coincides with 
onset of extra mixing was significantly overestimated by our 
models, especially in the low-metallicity regime. By consid¬ 
ering models of these two RGB mixing events we found that 
stellar codes better reproduce the structure of high metal¬ 
licity stars. To further examine this hypothesis we turned 
to photometric data of the MSTO and LF bump in a large 
sample of GCs. 

We directly compared the predicted MSTO and bump 
magnitudes as well as the parameter to the cor¬ 

responding values measured from photometry. We repro¬ 
duced two well known results: i) that 15 Gyr isochrones are 
required to match the MSTO in some clusters and ii) that 
the models match the bump luminosity at the highest metal- 
licities, but significantly overestimate it at low-metallicity. 
There is currently no reason to doubt the bolometric cor¬ 
rections at low-metallicity and there is no evidence to sug¬ 
gest that uncertainties in the metallicity scale are significant 
enough as to be responsible for the discrepancy. A system¬ 
atic underestimation of the distance modulus by a factor of 
1.03 would see 12 Gyr isochrones provide an upper limit to 
the ages of the GGs but would shift the models out of agree¬ 
ment with metal-rich clusters. A systematic error in the dis¬ 
tance modulus would not explain the fact that the models 
do not reproduce the observationally determined 
which is independent of distance and reddening. 

Changes to the stellar physics are the only means by 
which the theoretical value of AMv^^jP can be made con¬ 
sistent with observations. Our tests demonstrated that rea¬ 
sonable variations to the stellar microphysics do not signifi¬ 
cantly alter the stellar structure. Our investigations covered 
composition and a-element enhancement, opacity, convec¬ 
tive efficiency and the EOS choice on RGB evolution. Of 
these, updating MONSTAR to the most recent OPAL EOS 
yielded the greatest improvement (« 0.1 mag compared to 
the > 0.4 mag required). The inclusion of diffusive overshoot 
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in the models is sufficient to reproduce the bump magnitude 
and at intermediate and high metallicity. In or¬ 

der to directly match the bump magnitude at low-metallicity 
the efficiency of overshoot needs to be increased, so much so 
that agreement with is lost. The fact that our 

models with overshoot cannot simultaneously provide a limit 
(as a function of [M/H]) to the absolute magnitude of the 
bump and AMvb^i^jP is most likely due to predicted MSTO 
magnitudes that are too bright for 12 Gyr isochrones. 

The final test of the models focussed on their ability to 
reproduce the abundance changes associated with extra mix¬ 
ing. We determined that the thermohaline mechanism could 
not simultaneously account for the reduction of [C/Fe] and 
A(Li) as a function of RGB luminosity. As it stands, current 
observations imply that at low metallicity Li and carbon be¬ 
gin mixing a different luminosities. This result is indepen¬ 
dent of models and seemingly impossible. In the six GGs ex¬ 
amined, a free parameter between Ct = 120 — 300 provided 
a bounds to the Li data. This is in contrast to a value of 
Ct = 1000 cited throughout the literature to match the car¬ 
bon abundances in various stellar environments. Whilst this 
discrepancy is a challenge for thermohaline mixing, the con¬ 
straints may prove difficult for other mechanisms also, as Li 
is likely to be destroyed very efficiently under the conditions 
required to deplete C. Finally we suggest that the diffusive 
mixing algorithm used in stellar modelling may not be ade¬ 
quate to follow elements such as Li that are very sensitive to 
the local thermodynamic history. In order to quantitatively 
match observations numerical codes may require a more re¬ 
alistic treatment of mixing. 
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APPENDIX A: HEAVY ELEMENT MIXTURES 
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[a/Fe]^ 

=0.0 

[a/Fe] = 

=0.2 

[Q;/Fe] = 

=0.4 

[Ni/Nz] 

[i/Fe] 

[Ni/Nz] 

[i/Fe] 

[Ni/Nz] 

[i/Fe] 

c 

2.58x10-1 

0.0 

1.95x10-1 

0.0 

1.41x10-1 

0.0 

N 

6.49x10-2 

0.0 

4.90x10-2 

0.0 

3.53x10-2 

0.0 

O 

4.71x10-1 

0.0 

5.63x10-1 

0.2 

6.43x10-1 

0.4 

Ne 

8.17x10-2 

0.0 

6.17x10-2 

0.0 

4.45x10-2 

0.0 

Na 

1.67x10-3 

0.0 

1.26x10-3 

0.0 

9.09x10-1 

0.0 

Mg 

3.82x10-2 

0.0 

4.57x10-2 

0.2 

5.22x10-2 

0.4 

Al 

2.71x10-3 

0.0 

2.05x10-3 

0.0 

1.47x10-3 

0.0 

Si 

3.11x10-2 

0.0 

3.72x10-2 

0.2 

4.25x10-2 

0.4 

P 

2.47x10-'! 

0.0 

1.86x10-1 

0.0 

1.34x10-1 

0.0 

s 

1.27x10-2 

0.0 

1.52x10-2 

0.2 

1.73x10-2 

0.4 

Cl 

3.04x10-1 

0.0 

2.29x10-1 

0.0 

1.65x10-1 

0.0 

Ar 

2.41x10-3 

0.0 

1.82x10-3 

0.0 

1.31x10-3 

0.0 

K 

1.03x10-1 

0.0 

1.23x10-1 

0.2 

1.40x10-1 

0.4 

Ca 

2.10x10-3 

0.0 

2.52x10-3 

0.2 

2.87x10-3 

0.4 

Ti 

8.56x10-3 

0.0 

1.02x10-1 

0.2 

1.14x10-1 

0.4 

Cr 

4.20x10-1 

0.0 

3.17x10-1 

0.0 

2.28x10-1 

0.0 

Mn 

2.58x10-1 

0.0 

1.23x10-1 

-0.2 

5.59x10-3 

-0.4 

Fe 

3.04x10-2 

0.0 

2.29x10-2 

0.0 

1.65x10-2 

0.0 

Ni 

1.59x10-3 

0.0 

1.20x10-3 

0.0 

8.67x10-1 

0.0 


Table Al. Heavy element mixtures used in generating type 1 opacity tables. The abundances are based on lAsplund et ^ ll2009t) and 
their determination of the solar photospheric composition. 













